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Abstract 

A closed form of the multi-peakon solutions of the Camassa-Holm 
equation is found using a theorem of Stieltjes on continued fractions. 
An explicit formula is obtained for the scattering shifts. 
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Camassa and Holm @,0 introduced the strongly nonlinear equation 

13 11 
(1 - -D 2 )u t = -{u 2 ) x - -(u 2 x ) x - ~{uu xx ) x , (1) 

as a possible model for dispersive waves in shallow water and showed that it 
could formally be integrated using the spectral problem (D 2 + k 2 m — = 0, 
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where m = 2(1 — ^D 2 )u. We showed in a previous article that a Liouville 
transformation (|]J, equation (5.8)) maps this spectral problem to the string 
problem 

v"(y) + k 2 g(y)v(y) = 0, -1< y < 1; u(±l,z) = 0, (2) 

where m(x) = g'(y) sech 4 (x), and y = tanh rc. 

The Liouville transformation is independent of m, and m may assume 
both positive and negative values. The scattering data are invariant under 
the transformation, and consist of a discrete set of eigenvalues fc| and asso- 
ciated coupling constants Cj. The evolution of the data under the flow (TO) is 
kj = and dj = -2 Cj /k 2 ; gf. 

The n-peakon solution has the form 

n 

u(x,t) =^2pj(t)exp(-2\x-Xj(t)\), (3) 

3=1 

where the Pj, Xj evolve according to a completely integrable Hamiltonian 
system 0. A number of properties of the two-peakon and n-peakon solutions 
were obtained from an analysis of the dynamical system in 0. 

In this note we use a continued fraction expansion and a theorem of 
Stieltjes to give explicit algebraic formulas for the n-peakon solutions. 

The multi-peakon solution arises when g(y) dy is a sum of delta functions: 

n 

g(y) d v = Yl 9 ^^ y ~ dy > -i < yi < ■ ■ ■ <y n <i- (4) 

Equation (0) is then interpreted in the following sense: 

v " = °> Vj-i < y < Vj\ v '(yj+) - v '(yj-) = z m v (yj)> ( 5 ) 

for < j < n. We let y = — 1 and y n+ \ = 1, and we have taken z = —k 2 
as the spectral parameter. We denote the lengths of the subintervals by 
h = Vj+i ~ Vj> for j = 0, . . . , n. 

Let if and ip be the solutions of (^) that satisfy the boundary conditions 
(p(—l,z) = 0, <f'(—l,z) = 1, ip(l,z) = 0, and ip'(l,z) = —1. The eigenvalues 
{Xj} of (H) are the zeros of <^(1, z), which are simple. The coupling constants 
{cj} are characterized by the identities ip(y, Xj) = Cjip(y, Xj). Differentiating 
with respect to y and setting y = 1, we find that Cj = —ip'(l, Xj). 
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To construct ip we define qj = ip(yj,z), and pj = (p'(yj—,z). Thus Pj is 
the value of the derivative <p' in the interval (yj-i, Vj)- The pj and qj can be 
obtained recursively from (0) by q% — lo, p\ — 1, and 

9j - Qj-i = ij-i Pj , Pj ~ Pj-i = z gj-i qj-i. (6) 

It follows inductively that (p(l, z) is a polynomial of degree n. 

We assume for most of this note that the n point masses gj are positive, 
so that the eigenvalues are negative: A n < A„_i < . . . Ai < 0. By classical 
oscillation theory, the Cj alternate signs, with c\ > 0. 

The scattering data is encoded in the Weyl function 

w(z) = (f/(l,z)/(f(l,z). 

Since <p(y, 0) = 1 + y, <p(l, 0) = 2; and therefore (p(l, z) — 2 Hj=i(l ~~ z /^j)- 
Combining these results, we obtain 



^ i = ^ + E 7 ^. «* = ^ lid - V*> (7) 

j=l ■> k^j 



The residues cij of w(z)/z at z — Xj are positive, since both the coupling 
constants Cj and the successive products alternate signs. 

Following Krein H , we may write the Weyl function as a continued frac- 



tion 



(8) 



zg r , 



This follows by induction on @. We have pi = 1, q% = Iq, qi = qi + I1P2, 
p 2 =p 1 + zgil ] hence 

P2 P2 1 



<?2 k+P2h , 1 

1. + r 

+ 7- 
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Assuming (|3|) for n — 1 masses, and adding an additional mass at y n G 
(y n -!, 1), we have p n+1 = p n + zg n q n , q n+1 = q n + l n Pn+x] and 



Pn+l 



Pn+l 



Qn+1 Qn + InPn+l 



ln + 



zg n + 



Pn 



A classical result of Stieltjes || recovers the coefficients of the continued 
fraction (^) from the Laurent expansion of w(z)/z at infinity, obtained from 
(0) by expanding each (z — A 3 -) -1 : 



w(z) 



£ 

fc=0 



-l) fc A fc 



A, 



2 + 1^7=1 a j 



hu(-h k ^ k>i. 



k = 0; 



(9) 



Since the eigenvalues are negative and the aj are positive, each is positive. 
Stieltjes showed that such a Laurent series can be uniquely developed in a 
continued fraction 



>2k 



A 1 A 1 ' 



hz + 



)2k+l 



(A 



1\2 
k) 



A° A 
^k^k+i 



b 3 z + ... 

Here Aq = 1 = Aj and A°, A£, fc > 1, are the k x k minors of the Hankel 
matrix 



H 



(Aq a 1 a 2 ...\ 

A 1 A 2 A 3 .. 
A 2 A 3 A 4 .. 



V : 



/ 



whose upper left hand entries are, respectively, A Q and A\. 

Comparing this continued fraction with that for the Weyl function, we 
obtain 



A -A 

n-3 n-j+1 



9j 



n-j+l, 



^n-j+l^n-j 



(10) 
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Because the aj are positive, the A° are positive, k < n + 1, and the A\ 
are positive, k < n. With n eigenvalues and coupling coefficients, A° +2 
vanishes and the continued fraction terminates. The time dependence of lj 
and Qj under (p]) is determined explicitly from the evolution of the coupling 
coefficients, (0), and @. 

The multi-peakon solution is given by 

1 f°° 

u(x,t) = - / exp(— 2[rr — x'\)dM{x , £), 



where c?M = g(y)(dy/dx)dx is the pull-back of g dy under the transformation 
y = tanh x. Hence 

n 

dM(x,t) = / ]gj{t) sech 2 (xj(t)) 5(x — Xj(t)) dx. 

3=1 

The positions Xj(t) are obtained recursively from 

1, 1 + Vj S^j i 

2 1 - K 

Since 4 = 2 we obtain 

A 

2 ° A 



x^^logjL, Aj=$>, A+ = $>. (11) 



The n-peakon solution is thus given by (^) with = 2<7j'Aj At. 

These formulas imply that the asymptotic positions of the peaks are given 

by 



n _ j+1 ~t/A J - + ilog2a i (0) + 5^1og(^-l), t^-oo; (12) 

■~*/A i + ilog2a i (0) + ^ log(l-^) t^oo. (13) 

k>j 



From these formulas it can be determined that the phase shift in the peakon 
with speed 1/Aj as t — > ±oo is 

E'o g (i-^)-E^(^-0- (") 

k>j k<j 
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Moreover, the asymptotic height of this same peakon is —1/Xj. From this we 
may conclude that the n-peakon solution is asymptotically a sum of single 
peakons moving with speeds —1/Xj. 

The two-peakon solution has the form (|3|) with 

u\ X\ai + X\a 2 ai + a 2 , , 

= ~TT7\ 7~\ V' P2 ^ = —\ T~\ ; 15 ) 

AiA 2 (,Aiai + A2CI2J Aiai + X 2 a 2 

^iW = ^ lo g^T5 — A2 ! 2 aiQ2 ' x 2 {i) = \\og2(a l + a 2 ), (16) 

where dj{t) = aj(0)e 2 '/ Aj '. The explicit form for the relative positions and 
momenta x 2 — X\ and p 2 — Pi was given in [[|. 

We have concentrated here on the multi-peakon case, for which the masses 
Qj all have the same sign and the solutions do not develop singularities. 
The procedure is purely algebraic, however, so the formulas hold also for 
mixed positive and negative masses: the peakon-antipeakon case. In this 
case singularities can occur. 
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